The vortex pinning on the cylindrical defects and 
the electronic structure of the vortex core 

V. L. Kulinskii, D. Yu. Panchenko 

Department of Theoretical Physics, Odessa National University, Dvoryanskaya 2, 
65026 Odessa, Ukraine 

E-mail: kulinskij@onu.edu.ua, dpanchenko@onu.edu.ua 

Abstract. The model of the Abrikosov vortex pinning on a cylindrical defect is 
proposed. It is shown that in the limit k 1 the potential for the vortex core 
excitations can be treated in terms of the zero-range potentials method. Using the 
variational method the estimates for the energy of pinning, the pinning force and the 
density of critical current defect are obtained. 

Keywords: Abrikosov vortex, pinning, bound states, zeroth range potential. 
PACS numbers: 74.25. Qt, 74.25.Wx, 74.25.Ha 

Submitted to: J. Phys. A: Math. Theor. 
1. Introduction 

The pinning of Abrikosov vortices on the defects influences the static and dynamic 
properties of the superconductors. In particular the critical current density j c which 
determines the stability of the superconducting phase depends on the pinning. In order 
to get the superconductor with higher value of j c the defects of the crystal structure 
should be created since they are the pinning centers. If there are too much defects 
the superconductivity is destroyed 0, |[|. Evidently the pinning essentially depends on 
type of defects interacting with vortex. In this connection the pinning can be of the 
following types: single - vortex pinning on the microscopical (nanometer scale) inclusion 
of normal phase ||; superficial pinning on a film inhomogeneity [|]]; pinning on a lengthy 
flat defects of twin boundary type or antiphase boundary defect ||; pinning on a 
linear (quasionedimensional) defects, in particular on screw and edge dislocation [J7|, |8|, |9| ; 
The most effective pining occurs when the geometry of the defect coincides with the 
geometry of the vortex, i.e. for the columnar defect parallel to the axis of the vortex 
[TO . Due to progress in nanotechnology now it is possible to implant the material with 
the defects of controllable characteristics of atomic size scale. For the physical situation 
considered the most interesting is the perforating the material with the cylindrical pores 
of radius R via the bombardment by fast ions or highly energetic protons [|ll], 1(| . 
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The phenomenological approach to estimate the force (per unit length) and the 
corresponding density of critical current in the case when the radius of the cylinder 
R <C A was proposed in [12| and extended to the case of arbitrary large values of R in 
T3| . These approaches were based on the phenomenological Landau-Ginzburg theory 



and therefore use the assumptions about the order parameter A(r), which is the main 
characteristic of electronic structure of the vortex. The case of R < £ was considered 
within the quasiclassical approximation for the Bogolyubov-de Gennes (BdG) theory 
|14| . The characteristic scale of the spatial variation of A(r) is the correlation length 
£. Although in many works the specific spatial dependence of A(r) is not of principal 
interest the consideration of the transition to the limit x — > oo where two characteristic 
scales R and £ merge does depend on the assumption of the spatial distribution of the 
order parameter A. According to |15| the structure of A(r) can be even more complex 
and is characterized by additional scale £1 < £, which separates the regions at the point 
where the jump of the derivative of the order parameter occurs. The quantity £1 also 
determines the distance where the supercurrent density reaches its maximum |16, [17]]. 
So we treat the distance ^ as the characteristic length scale of the vortex core. In its 
turn the increase of the slope of A(r) in the region r < £1 is reflected in "Kramer- Pesch 
effect" [|18j. It diverges in the quantum limit T — » where kp 1 ~ £i <C ^bcs because of 
the oscillations of the Bogolyubov wave functions The shrinking of the core region 
leads to a reduction in the number of bound states [[L9| . 

Thus the structure of the vortex core is far from trivial and should be treated 
correctly, especially in the limit x — > oo due to point-like singularities in the spatial 
distribution of physical quantities. 

As the vortex core plays the role of the defect the appearance of the bound states 
is expected in analogy with the localized states near the defects in semiconductors 120 . 
There exists the branch of the localized excitations which can be called "magnetic". 
Another branch of bound states is possible due to the repulsive interaction with the 



defect in the vortex core 22 



The aim of this paper is to formulate the model for pinning which treats the 
appearance of the bound states in the limit x — > oo. This limiting case can be considered 
within the zero range potential method f23|. We also neglect the temperature effect 
as we work in the quantum limit [|TIJ. Here only lowest bound states are populated 
and the behavior of the slope of A(r) at the vortex center becomes very steep. The 
structure of the paper is as follows. In Section we consider the relation between 
the BdG Hamiltonian for the qusiparticle excitations and the Aharonov-Bohm (AB) 
Hamiltonian for the particle in the field of the localized magnetic flux. We show that 
these hamiltonians are equivalent for the low lying energy states localized near the 
vortex core where the contribution due to order parameter A can be neglected. We 
apply the treatment of the AB Hamiltonian developed in for the investigation of 
these bound states. Basing on the physical nature these states we differentiate between 
the boundary conditions corresponding to the AB effect and to the case of the Abrikosov 
vortex (Ab- vortex). In Section || two types of the bound states are considered and their 
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susceptibility to the magnetic flux is analyzed. The results are summarized in the 
concluding section. 

2. The Hamiltonian for the bounded states in the vicinity of localized 
vortex 

The pinning mechanism depends on the electronic structure of the Ab-vortex core, the 
type of a defect and the interaction between them. In simple situation such interaction 
is due to that between electrons of the vortex core and the defect. Specific form of 
the interaction potential is determined by the interaction between the excitations which 
form the normal Ab-vortex core and the condensate of the Cooper pairs. 

The Bogolyubov-de Gennes Hamiltonian (BdGH), which describes the energy levels 



for such excitations is |2T^ 

F = ^ j^(p-a z ~A-~^W^ -E F }+a x A(r) (1) 

where Ep is the Fermi energy, the vector potential A describes the applied magnetic 
field of order H Cl while the gradient term is for the magnetic field localized in the vortex 
|2§. 

Since 

A ~ H Cl r; eA/chVO ~ #£ 2 /$ ~ H/H C2 < 1 , 
then we can neglect A and rewrite the BdGH (|l|) in the form (see f21| ) 

2m I dr 2 r dr \ 2 / r 2 
where \x is the angular momentum. 

The order parameter A(r) is the spatially variable which has obvious asymptotic 
behavior: 

f if r -> 

A(r) = (3) 
I A if r — )• oo 

The limiting case x — > oo corresponds to the point-like vortex. As we mostly interested 
in localized states near the Ab-vortex core where A is suppressed the order parameter 
A-term in ([I]) can be replaced by the proper boundary conditions at r — > 0. The inner 
structure of the vortex is encoded into the parameters of the boundary conditions. In 
particular in such an approach the divergence of the slope dA/dr at r — > can be 
treated correctly via introducing the parameter ^ dA/dr\ r=Q . It controls the specific 
boundary conditions. 

Note that without A-potential the Hamiltonian ([1]) coincides with that used for 
the description of AB effect where the localization of the magnetic flux also takes place 
. The elementary excitation with e Aoo are localized near the Ab-vortex line and 



as was noted in |2T| play a major role in transport and relaxation phenomena at low 
temperatures. 
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It is clear that the presence of a defect influence the spectrum of these excitations. 
The applied magnetic field ~ H Cl and therefore the vector-potential is of order 
A ~ H Cl £. Therefore with account of 

(4) 

and jp- ~ - << 1 the applied magnetic field for the excitations of interest can be 



omitted: 



where 



H = H AB + a x A , (5) 



&AB = 0z {i ( p + a *- c A '") 2 ~ Ef ] (6) 

is the part of the Hamiltonian. It incorporates the action of the magnetic field due to 
the localized vortex 

Aeff ~ H c £ . (7) 

If we take into account point-like character of the magnetic flux distribution: 
$ 

Aeff(r) 



and use the standard separation of radial and angular variables we obtain the Aharonov- 
Bohm Hamiltonian (ABH). Its radial part corresponds to that of the BdGH with a = 1/2 
(0) in the form 

Hab = <t z —\ ~— 2 ---+//- aa z f — Ep \ (8) 
2m L dr z r dr r z ) 

where a = {<f^} is the fractional part of the magnetic flux $ in the vortex core and $o 
is the magnetic flux quantum. 

As is known |T^| in the quantum limit the size of the vortex core £i <C £ then the 



magnetic flux which localized in the vortex core is defined as $ = 7T^fH C2 . It is easy to 
show that a has the form: 
$ 1 

a = 

where 



mu Hc2 



is Ginzburg-Landau (GL) coherence length. 
Also 



2e# C2 
2 m c 



is the cyclotron frequency. 
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Based on the fact that in the quantum limit £1 



i 

kp 



HI we obtain 



a 



AE f 



(10) 



This suggests that energy of bound state Huh, 
section |3|) having regard to @ is proportional to (^j- 

bound states (H), then it is natural to measure the energy in units of (^-j E F . 
If one takes into account that 



more detail it will be discussed in 

2 

E F . Thus the equation for the 
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then 



where 



1 

a 2 



dA 



dr 



r=0 



TT 



7T 



4/c, 



dr 



is the BCS coherence length, A the gap at T = and vp is the Fermi wave velocity. 
According to the work |19| for high-temperature superconductors in which 1 < kp^ < 5 
in the quantum limit the value £i/£o lies in the interval 0.3 -j- 0.6. Therefore a <C \ that 
leads to the fact that the magnetic flux $ in the vortex core is very small compared 
with $ - So for example in the quantum limit for the YBCO (kp£o ~ 4 with T c = 90 K 
|27|) 5 fi/fo ~ 0.3 that corresponds to a ~ 0.05. 

The Hamiltonian (|j) can be treated by the zero-range potential method [^| 
by specifying the correct boundary conditions which provide self-conjugacy of the 
Hamiltonian (^). In |23J] all possible boundary conditions for the Hamiltonian (||) as 
well as its spectrum were obtained. The general result is that depending on the type 
of the boundary conditions there can be no more than two bound states in the singular 
point where the magnetic flux is localized. Although AB effect is interesting in itself 
it can not provide clear physical interpretation for all types of the boundary conditions 
because the only output is the interference shift. 

We show that the existence of the bounded states of the elementary excitations in 
the the Ab-vortex and the pinning defects can be described within such approach. We 
search for the interpretation of these nonstandard boundary conditions found in |24| in 
terms of the vortex pinning because of the close relation between these problems stated 
above in the considered limit. 



3. Two types of bound states near the vortex core 



In [ 21 1 the existence of the low lying branch of bounded states which is emerged from 
the vacuum state of Cooper pairs was proved. Because of the magnetic nature of the 
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vortex the low lying bound states can be thought of as Landau levels in the effective 
magnetic field H e ff ~ H C2 with the spectrum: 

E^ = hu Heff fx, (12) 

where 



1 dA 

WH etf = WL + — 



(13) 

o 



kp dr 

In extreme type-II superconductors the first term can be omitted in comparison with 
the second term which is much greater P2] ]. Note that the result ( ]T3] ) does not depend 
on the specific form of A(r) but only on the limiting value of the radial derivative. This 
is quite natural as the bound states are strongly localized. Therefore in order to get the 
ground state for this branch we can omit the A-potential term in ([!]). Then the proper 
boundary conditions which correspond to the localization of the magnetic flux and in 
particular the singular behavior of 4= Hamiltonian (§) should be used. 

According to the theory of self-conjugate extensions for the Hamiltonian @ [[24 



there is 4-parametric set of boundary conditions for each value of a. They are 
represented by the unitary matrix: 



l UJ 



(!('■■• y/l — q^e 



U = e , 

'1 - q 2 e~ lV2 qe~ llpl ) . (14) 
0<g<l, z = l,2 

Note that for the Hamiltonian (||) a stands for the value of the localized flux. In the 
case of the Ab-vortex this means that this is the fraction of the flux quantum in the 
vortex core. 

Regardless of the fact that Hamiltonians (§) and (|6]) coincide, the boundary 
conditions which determine their spectrum are related to different physical situations. 
Depending on the type of the boundary condition there can be one, two or none bound 
states J24|. The case of the absence of the bound states corresponds to the AB effect 
itself (q — 1, uj = 0, a = tt, a G (0, 1)). 



The result of |2T| allows one to state that in the case of localized Ab-vortex at least 
one bound state should exist. The bound state which may occur in the AB effect due 
to the point-like interaction with the material of the solenoid is different in physical 
nature from that which appear as the bound state of the quasi-particle excitation near 
the vortex core. These states should be different with respect to their sensitivity to the 
change of the magnetic flux a. We consider this case first. 



3.1. Bound states in the vicinity of the vortex core 

The existence of the only bound state can be related to 1) the AB plus the defect which 
leads to the localization by its point-like potential (g < 1, u = 0, a = tt, a G (0, 1)); 
2) the localized Ab-vortex (0 < q < 1, oj = 0, a — tt/2, a G (0,1/2)) where the 
lowest bound state is due to localization of the excitation near the core surrounded 
by the Cooper pairs. It is possible to distinguish between these situations considering 
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the dependence of the bound state on the change in magnetic flux ratio a. In case of 
AB bound state which is caused mainly by the defect the energy of this state does not 
depend essentially on a. In contrast to this the bound state for Ab-vortex depends 
on the flux ratio much stronger. To confirm this statement and to connect the case of 
the only bound state for the Hamiltonian (0) with the result ( pj|) let us consider the 
derivative which is the measure of the "magnetic sensitivity" of the corresponding 
energy level. Note, that both parameters a and /x enter the Hamiltonian by the same 
way. Therefore the value |^ for the Hamiltonian (H) can be compared in the quantum 
limit with the value 



dEjE F 



hw L /E F + ( ^ 



for the BdGH. Equation ([15]) also can be written as 
defmition of a in (|9[) we can obtain that a 



E F 



(15) 



2 and using the 



j , . . This is consistent with (|10|) 
To obtain the dependence q(a, rj) we use the condition 



£i ( oi, q 



V 



(16) 



to fix the energy gap. Here 

2 

v 



hw L /E F + ( ^° 



in accordance with ([L5|) . In accordance with said above the value of the derivative ^ 
as: 



da I Aq 



dei 
da 



A 



dei 
dq 



Ao,a(g) 



dq 
da 



(17) 



The dependence of ^ | Aq and a on the parameter q is shown in figure [1] and figure ^| 
are respectively. Note that the value |p- increases linearly with the parameter q and the 



value a is also increases linearly at low q. According to (|TT|)-(|T3D the same behavior ^ i 
and aon^ takes place. It allows to state that the parameter q controls the limiting 



dA 
1 dr 



at x — > oo. 



slope parameter £ 

It should be noted that the range of values q(a, rf) depends essentially on the 
parameter rj. There is the minimum value of rj, when q runs over the interval (0,1) 
if < a < 1/2. This minimum value is ^( mm ) = 2 ~ 2.11 , where W is Lambert 

function. This result is in correspondence with the following facts. The first one is that 
the existence of half-quantum vortices bounds the value of a by 1/2. The second is that 
in the quantum limit the value rj is of order of unit and can not be too small. In the 
opposite case f^- ~ 1 j when rj <C 1 there is no dependence between q and a, that leads 
0. 

The choice of the parameter a influence the bifurcation between one and two bound 
states, which occurs with the variation of the flux parameter a. Thus in the case of 
one bound state, the value of parameter a is associated with the maximum value of 



aa I Aq 
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the parameter a for this bound state. This corresponds to the maximum value of the 
magnetic flux in the Ab-vortex core. The total flux is distributed between the core of 
size £1 and the outer region where the magnetic field decay exponentially. As was shown 
28| in epitaxial films of Bi 2 Sr 2 CaCu20$ + s, a t the meeting point of a tricrystal substrate 



of SrTiO-j appears half-quantum vortex (HQV). This effect is associated with d-wave 
nature of the superconductor and together with photoemission results, proves that the 
in-plane order parameter for this high-Tc cuprate superconductor closely follows d x 2_ y 2 
symmetry. The existence of the half-quantum vortices allows to suggests that within the 
vortex core should be no more than half of magnetic flux quantum. Such a restriction on 
the fraction of the magnetic flux in the vortex core corresponds to the choice of a — |. 

3.2. Bound states generated by the pinning center 

The low lying bound states considered above appear due to the interaction with 
surrounding Cooper pairs which prohibit the propagation of the excitations beyond 
the core. If the Ab-vortex is pinned by the defect one can expect the appearance of the 
new bound state. From the physical point of view such bound states may appear as the 
result of the resonant scattering from the defect from one side and the Cooper pairs in 



the bulk of the superconductor from the other side. In |L4[ it was shown that indeed the 
additional branch of the bound states appears. If the results in [14j] rewritten in terms 
of /j,, then we can show that this branch is less "sensitive" to magnetic field, which by its 
nature is different from the bound states discussed above. This difference is concluded in 
the behavior of the derivative d E ^ Ef , which is less than zero. To compare these results 
with our consider the case when u > \ with fixed a = n/2 which leads to appearance 
of two roots. The presence of the second bound state is connected with the appearance 
of additional branch associated with the resonant scattering from the insulating defect. 
To confirm this statement let us consider the derivative ^ I . which is the measure of 

da IAo 

the "magnetic sensitivity" of the corresponding energy level. It is easy to see that this 
derivative is not zero in contrast to the previous case. The dependence of ^| Ao on 



the parameter a is shown in figure [| where shows that 
parameter a that is consistent with behavior - ■ 

— ^. _j i i i i i i i i i i i i i i \ 



^r 1 decreases linearly with the 



da 01 

- 0.01 : 

- 0.02 \ 

- 0.03 ': 

- 0.04 \ 

- 0.05 

- 0.06 \ 

- 0.07 \ 

Figure 3. The dependence of 



0.2 



0.3 



da Ia o 



on a. 



Conclusion 

In the paper it is shown that the consideration of the low lying bound states localized 



in the vicinity of the vortex core can be investigated with the help of the formalism p4 
developed for ABH. We show that the nonstandard boundary conditions constructed in 
|24| can be interpreted in terms of the localization of the excitations in the vortex core. 



It is shown that at least one bound state exists in the vortex core even in the quantum 
limit when £i ~ kp . Within the proposed approach the case of the divergent slope of 
the order parameter can be treated correctly Another possibility for the localization 
can be realized for the pinned vortex due to the resonant scattering of the excitation 
between the defect and the surrounding Cooper pairs. 

Appendix A. Characteristics of pinning 

Here we give the variational estimates for the characteristics of pinning. The result 
(O) allows to consider the lowest bound states in the core as the Landau levels in the 



effective magnetic field of order H C2 corresponding to the localized flux. The effect of 
pinning can be considered as the additional localization due to the defect. This gives 
the ground to choose the trial wave function of the ground state as the superposition of 
bound states for these potentials correspondingly. 

We choose the trial wave function of the ground state (m = 0) in the form: 

* = ^ Land Sin /3 + ^ bound COS f3 (A.l) 
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where 

^ Land 

and a cyclotron radius 



exp( 



p ' 

'4a 2 - 



bound 



Up) 



(A.2) 



2uH r 



are the corresponding wave functions representing the ground states for corresponding 
interaction. Pinning energy per unit length calculated based on a variational method 
for YBaCuO (£ ~ 12 A, A ~ 1000 A p9[|). The result for the energy as the function 
of the variational parameter is shown in figure [AT. 



"pin 



1V , cm 




Figure Al. The dependence of e P i n on the variational parameter /3. 



The corresponding minimum of the energy is: 

e pin = -1.45 • 10- 6 Erg/cm (A.3) 
The force of the pinning per unit length is calculated as: 

(A.4) 



Spi n 
Jpin= ~QR~ 

ie corresponding 

values 



R=i 



and the corresponding critical supercurrent j = jrfpin at the used parameters have 



f pm *l8^ } ^9-10 8 A 
cm cm z 

Using such a simple model the temperature dependence of e pin and j can be 
obtained. Indeed, according to the equations of (0); ( |A-.4j ) and the definition of j 
for the quantities e P i n and j the following relations are valid: 

3 

£pin H C2 ; j c oc Hc 2 ■ 
Next, we use the fact that H C2 ~ r, where r = 1 — (T/T c ) 2 and then we easily find that 

3 

(A.6) 

Note that such temperature dependence of j coincides with the result of [1211 . The 



calculations of e P i n performed for YBaCuO in give e P i n ~ —10 6 Erg /cm. This is 



consistent with the estimates (|A.3|) obtained above. 
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